Large nonlinear Kerr effect in graphene 
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Under strong laser illumination, few-layer graphene exhibits both a transmittance increase due to saturable 
absorption and a nonlinear phase shift. Here, we unambiguously distinguish these two nonlinear optical effects 
and identify both real and imaginary parts of the complex nonlinear refractive index of graphene. We show that 
graphene possesses a giant nonlinear refractive index n2 — 10~^ cm^W"-*^, almost nine orders of magnitude 
larger than bulk dielectrics. We find that the nonlinear refractive index decreases with increasing excitation 
flux but slower than the absorption. This suggests that graphene may be a very promising nonlinear medium, 
paving the way for graphene-based nonlinear photonics. © 2013 Optical Society of America 
OCIS codes: 160.4330, 160.4236, 190.3270, 190.7110. 



Graphene, a single sheet of carbon in a hexagonal lattice, 
exhibits many interesting electrical and optical proper- 
ties [1,2] which arise due to its particular bandgap struc- 
ture E±{p) = V^IpI, where the sign corresponds to elec- 
tron (respectively hole) band, p is the quasi- momentum, 
and V::^ 10^ ms~^ is the Fermi velocity. In particular, a 
single graphene layer absorbs a fraction 7ra ::^ 2.3% of 
the incident light across the infrared and visible range, 
where a is the fine structure constant [3]. 

Graphene also has a broadband ultra-fast saturable 
absorption [4-7] , due to valence band depletion and con- 
ductance band filling. This saturable absorption has been 
used, amongst other applications, for laser mode lock- 
ing [4,7]. 

The massless bandgap structure of graphene has led to 
the prediction of other strong nonlinear optical effects, 
including frequency multiplication [8, 9] and frequency 
mixing [10]. Recently, broadband four- wave mixing in 
few-layer graphene samples has been reported [11]. This 
allowed the determination of the absolute value of the 
third order susceptibility of a single graphene layer 
\Xgr I — 1.5-10"^ esu, approximately eight orders of mag- 
nitude larger than in bulk dielectrics. Four- wave mixing 
in graphene deposited on optical ferrules has also been 
observed [12]. 

To our knowledge, there has so far not been any 
measurement that unambiguously discriminates between 
the real and imaginary part of the third order susceptibil- 
ity x^^^ of graphene. We report herein the experimental 
measurement of the nonlinear optical refractive index n2 
of loosely stacked few-layer graphene using Z-scan tech- 
nique. As demonstrated in [13], the knowledge of this 
coefficient together with the nonlinear absorption is suf- 
ficient to calculate the complex value of x*^^^- 



Samples of graphene were prepared as follows: a few- 
layer graphene thin film was grown on 25 jiia thick Cu 
foils by chemical vapor deposition (CVD) [14]. The cop- 
per was etched with ammonium persulfate (0.1 M). After 
thorough rinsing in distilled water, the fioating graphene 
film was fished by quartz substrate and dried gently in 
nitrogen gas. The optical image of a typical graphene 
sample is shown in Fig. 1 (a). The non-uniform color 
contrast of the optical micrograph indicates that the film 
has a variable thickness, i.e.., the dark contrast zones cor- 
respond to thicker graphene domains whereas the faint 
contrast zones correspond to thinner graphene domains. 
This is further confirmed by the Raman map of the G 
band, as shown in Fig. 1 (b), in which bright zones have 
larger G band intensity, indicating more graphene lay- 
ers. Raman spectra given in Fig. 1 (c) are captured from 
marked spots in Fig. 1 (a), (b), with zone A bilayer and 
zone B 4-layer graphene. The thickness can be unam- 
biguously identified by optical contrast (not shown) on 
optical transmission measurements. Different graphene 
samples, with respectively 1, 2, 3, 4, and 6 layers were 
used in the experiments reported below. The fabrication 
process produces loosely stacked graphene layers wherein 
the optical properties are expected to arise from the cu- 
mulative contribution of each layer. 

We measured the real and imaginary parts of the com- 
plex nonlinear refractive index by using the Z-scan tech- 
nique [15]. Our experimental setup has been used previ- 
ously in [16], but has been adapted to the small thick- 
ness of the sample. The sample is subjected to picosec- 
ond pulses emitted from a Pritel picosecond laser with 
center wavelength 1550 nm, pulse duration 3.8 ps (full 
width at half maximum), 3-dB spectral width 1.2 nm 
and repetition rate 10 MHz. These pulses are ampli- 
fied through a Keopsys erbium-doped fiber amplifier. Af- 
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Fig. 1. Characterizationofgraphene samples, (a) Optical 
image of graphene sample on silica substrate. Scale bar: 
10 /im. (b) Raman map of the G (1560 to 1620 cm~^) 
band of the area in (a) (WITec alpha, Laser wavelength: 
532 nm, spot size: ~ 500 nm, 100 x objective). CCD cts.: 
charge-coupled device counts. Scale bar: 10 /im. (c) Ra- 
man spectra from the marked spots: A (B) indicates 
the area with fainter (darker) contrast given by thinner 
(thicker) graphene. 



Fig. 2. Z-scan traces for samples with three-layer 
graphene taken at average power of 3 mW, corresponding 
to a peak power at focus of 0.56 GWcm~^. (a) Near field 
(open aperture), (b) Far field (closed aperture). Upon 
dividing by the near field curve one obtains the data of 
panel (c) which exhibits the typical shape of a Z-scan 
curve with positive nonlinear phase shift having an on- 
axis value of A<l> = 0.147 rad. The dotted curve was 
recorded in cw at the same wavelength and mean power. 



ter out-coupling from the fiber, shaping and routing the 
beam through the focusing objective, the average power 
of the pulse train can be adjusted between and 4.5 mW. 
From the laser parameters, the peak pulse power can 
therefore be adjusted from to 118 W. The beam is 
focused using a 20x microscope objective, generating a 
beam waist of 3 /im, corresponding to illuminations rang- 
ing up to a peak intensity / of 0.84 GW cm~^. The beam 
waist was measured by the beam profiler and further con- 
firmed by z-scan parameter fitting. The graphene sample 
and its quartz substrate are oriented perpendicularly to 
the beam axis and translated along the axis through the 
focus with a linear motorized stage. A portion of the am- 
plifier output is picked off and measured to monitor the 
optical power level continuously. 

Measurements are performed in two regimes, an 
open-aperture regime wherein all the light transmitted 
through the sample is collected on a photo-detector and 
a closed-aperture regime where only an on-axis portion 
of the diffracted beam is collected. The open-aperture 
regime enables the characterization of the intensity- 
dependent absorption. A typical trace when the sam- 
ple is translated through the beam focus is shown in 
Fig. 2 (a). The peak in transmission when the sam- 
ple passes through the focus is characteristic of sat- 
urable absorption. The sample's linear absorption al- 
most completely vanishes for illuminations greater than 
0.6 GWcm~^. Note that the wings visible in Fig. 2 (a) 
present an unusual slope which may arise from the sam- 
ple's inhomogeneity. Its presence does not affect our con- 
clusions. 

A typical closed aperture measurement is shown in 
Fig. 2 (b). In this trace the effect of the nonlinear phase 
is of the same order of magnitude as the effect of sat- 



urable absorption. In order to isolate the former, we di- 
vide the curve in Fig. 2 (b) by the curve in Fig. 2 (a), 
thereby obtaining the curve in Fig. 2 (c). The latter 
has the typical shape of a Z-scan trace. The up-down 
curve implies a small positive on-axis phase shift A<l>. 
For small A<l>, one can fit the curve with the function 
T{x) = 1 + (i+^?)f9+a;2) where x = -z/zr is the nor- 
malized distance from the focus and zr is the Rayleigh 
length. The above method for deriving A^ from the 
measured data is applicable when the distortion of the 
phase front is small and the closed aperture is sufficiently 
narrow [15] . We also performed a full calculation of local 
nonhnear effects across the graphene layer and its sub- 
sequent propagation, and then fitted Fig. 2 (b) directly. 
Both approaches provide comparable results. Additional 
measurements carried out in cw regime with similar av- 
erage power as in the pulsed regime {c^ 1 mW) show that 
our Z-scan measurements are not affected by cumulative 
thermal effects (see Fig. 2(c)). This is presumably due to 
the very high thermal diffusion coefficient of graphene. 

Z-scan measurements under variable optical power 
were performed. Closed-aperture measurement results 
show that the transmittance difference AT, which is 
the difference between the peak and the background in 
Fig. 2 (a), has a power dependence characteristic of sat- 
urable absorption. We fitted the transmittance difference 
by AT = ATo- ATo/(l + ///sat), see Fig. 3 (a), yielding 
the estimate /gat = 74 MWcm~^, which is of the same 
order of magnitude as the value reported in [17], but 
an order of magnitude smaller than the values reported 
in [5] . From this value of /gat , we can calculate the non- 
linear absorption coefficient P that appears in [13]. The 
maximum transmittance difference for 3-layer graphene 
is ATo = 5.1%, indicating almost complete saturation of 
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Fig. 3. (a) Relation between transmittance difference 
(AT) and input power for three layer graphene. (b) De- 
pendence of A<l> (right axis) and n2 (left axis) on peak 
intensity for three layer graphene. (c) Relation between 
A<l> and number of layers, for different peak intensity. 



the absorption. Measurement of the maximal transmit- 
tance difference ATq shows that it scales linearly with 
the number of layers (not shown). 

The nonlinear phase A<l> is plotted on Fig. 3 (b), with 
respect to the input power. The Kerr refractive index n2 
can be deduced from the slope of this curve at low inten- 
sities, using 722 = A^/{koLI) where ko = 27r/A and L 
is the sample thickness, assumed equal to 1 nm. A value 
of n2 — 10~^cm^W~^ is obtained, which is approxi- 
mately 10^ times larger than that of bulk dielectrics. 
When the intensity is increased, starting from / > Igat ^ 
O.lGWcm"^, the change in the nonlinear phase A$ 
saturates. In this high intensity regime, the exact mod- 
elling of the nonlinear response requires to take into ac- 
count x^^^ and higher order odd terms, such as x^^\ • • • 
However, from the experimental data, an effective non- 
linear index can be defined by n2(/) = A^/{koLI). 
One sees that nj decreases with increasing /, and then 



reaches a constant value n2^Q- 10~^cm^W~^ for 
/ > 0.6 GW cm~^. We also measured the nonlinear phase 
as a function of the number of graphene layers and found 
that at fixed powerit scales linearly with the number of 
layers, see Fig. 3 (c). The real part of the complex nonlin- 
ear refractive index of loosely stacked graphene is found 
to be independent of the number of layers. 

Note that the value of x*^^^ based on four- wave mixing 
experiments reported in [11] corresponds to an equiva- 
lent 722 of 1.5-10"^ cm^ W~^, which is forty times smaller 
than our value. We attribute this discrepancy to the fact 
that the third order susceptibility of graphene proba- 
bly has multiple origins. The parametric process of four 
wave mixing measured in [11] arises from the coherent 
electronic response and the massless bandgap structure 
of graphene, as explained in Refs. 8-10. In addition there 
is a strong non-parametric process that gives rise to sat- 
urable absorption, and which presumably also gives rise 
to a strong non-parametric contribution to n2 . For a gen- 
eral discussion of these processes, see [18, p. 13] The Z- 



scan measurements reported here measure the sum of all 
contributions to n2. 

In conclusion, the characterization of the nonlinear 
optical properties of loosely stacked graphene was per- 
formed on a Z-scan setup. Measurements of the saturable 
absorption are in agreeement with previously reported 
values. Measurements of the nonlinear phase yield a non- 
linear coefficient n2 — 10~^cm^W~^, where n2 is de- 
fined in [13]. The measured value of the effective non- 
linear refractive index 712 decreases with increasing in- 
put power. Interestingly it remains large when the ab- 
sorption is saturated. This suggests that at high power 
(/ > 0.6GWcm~^) the figure of merit of graphene for 
nonlinear optics applications may be highly favorable 
(high nonlinear phase, low absorption per layer). 
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